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Abstract 

We investigate the propagation of electromagnetic waves in resistive e^ pair plasmas using a one- 
fluid theory derived from the relativistic two-fluid equations. When the resistivity normalized by the 
electron/positron inertia variable exceeds a critical value, the dispersion relation for electromagnetic 
waves shows that the group velocity is larger than the light speed in vacuum. However, in such a 
case, it also is found that the plasma parameter is less than unity: that is, the electron-positron 
pair medium no longer can be treated as plasma. Thus the simple two-fluid approximation is 
invalid. This confirms that superluminal propagation of electromagnetic wave is forbidden in a 

plasma a conclusion consistent with the relativistic principle of causality. As an alternative, 

we propose a new set of equations for "causal relativistic magnetohydrodynamics", which both 
have non-zero resistivity and yet are consistent with the causality principle. 

PACS numbers: 52.27.Ny, 52.30. Cv, 52.30.Ex, 52.35.Hr 
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I. INTRODUCTION 



Relativistic magnetohydrodynamic (RMHD)numerical simulations have been performed 



lOl, 



11 



12l. |13H. These numerical 



by a number of groups recently 
simulations revealed many important, interesting features of relativistic plasmas, especially 
around rotating black holes in active galactic nuclei (AGNs), microquasars, and gamma-ray 
bursts. Regarding energy extraction from a rotating black hole, the magnetohydrodynamic 
(MHD) Penrose process has been confirmed [4, 5], and long-term simulations of relativistic 



jet formation around a rotating black hole have been performed [13|. All of these RMHD 
simulations were restricted by the ideal MHD condition, where electric resistivity is zero. 
In spite of recent significant advancements in ideal RMHD simulations, one with finite 
resistivity (resistive RMHD) have not been performed seriously except in a few cases (e.g. 
Q, I5I). This is reasonable because there has been concern that the inclusion of finite 
resistivity in the RMHD equations destroys their causality. In fact, the group velocity of 
electromagnetic waves derived mathematically from the resistive RMHD equations is larger 
than the light speed in vacuum. This raises the possibility of superluminal communication, 
which is contradictory to the relativistic principle of causality. The main purpose of this 
paper is to clarify and rectify this problem. 

To fix this problem, we must reconsider the resistive RMHD equations. Such a task was 

n 

first performed by 161] using the Vlasov-Boltzmann equation for a pulsar magnetosphere. It 
yielded a relativistic version of the generalized Ohm's law and a new condition for the validity 
of the MHD approximation for a pulsar magnetosphere (where the Lorentz factor is much 
larger than unity). A more generalized treatment, which included annihilation of electrons 
and positrons, radiation, Compton scattering, and pair photoproduction was formulated by 



171] and 18|]. Reconsideration of ideal MHD in a neutral cold plasma based on two-fiuid 



approximation was presented by [19| , who investigated the conditions under which the MHD 
approximation breaks down. For investigation of black hole magnetospheres, formulated 
the general relativistic version of the two-fiuid approximation in the Kerr metric. An even 
more generalized version in a time- varying space-time was derived by 2l| from the general 
relativistic Vlasov-Boltzmann equation. 

In this present paper, we derive the one-fiuid equations of an electron-positron (pair) 
plasma based on the two-fiuid equations with a new definition of variable averaging for the 
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two fluids (section [IT])- In section UTTt we derive the dispersion relation for electromagnetic 
waves in uniform, unmagnetized and magnetized pair plasmas. We then examine the sit- 
uation where the group velocity of electromagnetic waves in the resistive plasma is larger 
than the light speed in vacuum and show that this condition cannot be realized in a plasma 
whose plasma parameter is larger than unity. In section \IV\ we propose a simple set of 
resistive RMHD equations, which are consistent with the principle of causality. In section 

rV] we discuss phenomena with respect to the superluminal propagation of wave packets 

phenomena that cannot be avoided when the RMHD equations are acausal. Finally, our 
summary and discussion are presented in section IVll 



II. RELATIVISTIC TWO FLUID MODEL OF PAIR PLASMA 



To provide a solid base for resistive RMHD, we begin with a relativistic two-fluid model 
of a pair plasma in the Minkowski space-time x^, x^, x^) = {t,x,y,z), where the line 
element is given by ds'^ = —{dx^Y + (dx^)'^ + (c/x^)^ + (rfx^)^ = rj^^dx^dx^ . Throughout this 
paper (except for one paragraph in section UTTl) . we use units in which the light speed, the 
dielectric constant, and the magnetic permeability in vacuum all are unity: c = 1, eo = 1, 



/io = L 



[e.g., 
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he relativistic equations of the electron and positron fluids are given as follows 



23|): 



^^.{n^U^) = 0, (1) 
d^{h^U^U-^) = ± en^r^^'^U^F,, ± i^^ (2) 

d, *F^'' = 0, (3) 

d^F^" = r, (4) 

where a variable with subscript, plus (+) or minus (-), is that of the positron and electron 
fluid, respectively, n± is the proper particle number density, f/^ is the four-velocity, e is 
the electric charge of positron, p± is the proper pressure, h± is the relativistic enthalpy 
density, F^,^ is the electromagnetic fleld tensor, *F^'-' is the dual tensor density of F^^, 
is the frictional four-force density between the electron and positron fluids, and is the 
four-current density. We often will write a set of the spacial components of the four-vector 
using a bold itahc font, e.g., U± = {U]_,U^,U^), J = {J^,J^,J^). Here we assume that 
the electron/positron fluids are heated only by Ohmic heating and neglect pair creation and 



annihilation. We also neglect radiation and quantum effects. 

We further define the Lorentz factor 7± = f/±, the three- velocity VJ^ = f/±/7±, the electric 
field Ei = the magnetic flux density Bi = ^ Ylijk^ijk^^^ Levi-Civita tensor), 

and the electric charge density pc = ■ Here, the alphabetic index (i, j, k) runs from 1 to 3. 
Using the above relativistic equations we obtain the vector form of the relativistic 

two-fluid equations, 

^(7±n±) + V-(n±C/±) = 0, (5) 

^ {h±U±) + V ■ {h±U±U±)=-Vp± ± e-f±n±{E + V±x B)±R, (6) 



dt 

^(7i/^± - P±) + V • (7±/i±C/±) = ±en^U ■E±W, (7) 

V ■ = pe = e(7+n+ - 7_n_), (8) 

V • B = 0, (9) 

^ ^ -V X (10) 

— + J = VxB. (11) 

The frictional four-force density between electrons and positrons is (equation (1A9P in Ap- 
pendix , 

-R^ = [n^'y^n^UI — n^j,n^U_] H 7^ —it , (12) 

n ^+7+ + n_ji_ 

where Uee is the electron-positron Coulomb collision frequency, m is the mass of an elec- 
tron/positron particle, the variables with primes are physical quantity observed in the center- 
of-mass frame of the two fluids, and R^' is the energy gain rate of the positron fluid due 
to the friction with electron fluid in the center-of-mass frame. Note that the variables ob- 
served in the center-of-mass frame are proper variables (see Appendix . If we assume that 
relative velocity of the positron and electron fluids is much smaller than thermal velocity 
of the fluids, the collision frequency i^ee is proportional to the relative velocity of the two 
fluids and its coefficient depends only on temperature and density of the two fluids. On the 
other hand, if the relative velocity of the positron and electron fluids is relativistic, then the 
coefficient also depends on their relative velocity. Note also that when the relative velocity 
of the electron and positron fluids is nonrelativistic, 'j'^ ^ 1, while in the case of relativistic 
relative velocity, 7^|_ > 1. Through this paper, we usually assume that the relative velocity 
of the two fluids is smaller than the sound velocity of the plasma (which is nonrelativistic). 



To derive the one-fluid equations of a pair plasma, we define average and difference 
variables as follows: 





n+ + n„ 
2 ' 




P = 


2mn, 


(14) 


= 




(15) 


2n 


J^" = 


e{n+U^-n_U^), 


(16) 


h = 




(17) 


Ah = 


\nz^ n_) 


(18) 



From the relativistic two- fluid model of the pair plasma (!5l)-( [T2l) . we then can obtain the 
one-fluid equations of the pair plasma, 



^(7P) + V ■ (pC/) = 0, 
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V-E = Pe, (23) 

V ■ S = 0, (24) 

f ^ -V X E. (25) 

— + J = V X B, (26) 

where rj = mvee/{ne^) is the electric resistivity and is the equipartition factor of the 
thermal energy due to the friction between the electron and positron fluids given in Appendix 

m 



Equation (12T!) corresponds to the generalized Ohm's law. The classical formulation of the 
resistivity r/ of a non-magnetized nonrelativistic plasma (due to Coulomb collision) is used, 
and the formulation for an electron-proton plasma is confirmed as an appropriate expression 
for resistivity of weakly magnetized nonrelativistic plasma in laboratory experiments (e.g., 



in a "Tokamak" thermonuclear fusion device 
pair plasma. 



2J|). Note that the Hall effect disappears in a 



III. DISPERSION RELATION FOR ELECTROMAGNETIC WAVES IN RESIS- 
TIVE PLASMA 

In this section, we derive the dispersion relation for electromagnetic waves in a pair 
plasmas using a linear analysis of equations (fT9|) - (l26|) . First, we assume that the background 
plasma is at rest, uniform, and non-magnetized: p = Pq {n = Uq), V = 0, p = po {p+ = p_ = 
Po/2), h = ho, B = E = 0. Perturbations due to the electromagnetic waves are so small 
that the plasma motion is non-relativistic. Then we have the following linearized equations 
with respect to the perturbations, pi = p — po {rii = n — no), Vi = V, pi = p — po, Bi = B, 
El = E, Ji = J , Pel = Pe, hi = h — ho, 

(27) 

(28) 

(29) 

(30) 
(31) 
(32) 
(33) 

(34) 



d_ 
di 



d 


+ V-{poVi) 


= 0, 


dt 


= -Vpi, 




El - 7] J I 


dJi 


{hi - pi 


)+V-{hVi) 


= 0, 




V-Ei 


= Pel, 




VBi 


= 0, 




OB I 

dt 


= -V X ^; 




dt ^-^^ 


= Vx Bi, 



where k = hQ/{2nQe)'^, and we assume p\_ p\. These linearized equations do not depend 
on the equipartition fraction of frictionally thermalized energy, ^ (0 < < 1). When we 
consider transverse modes of the linearized equations, fc • Vi = 0, the dispersion relation for 
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the electromagnetic waves can be written as 

where k is the wave number vector, uj the frequency of the electromagnetic wave, and 
t'ee — v/i'^'^) = ("^'^o/^o)'^cc- Here we also obtain pi = 0, pi = 0, pci = 0. 

When we normalize the variables as a) = uj/{2v'^^ and k = fc/(2z/gg), the dispersion 
relation (l35l) becomes 

H{u^ - P) (1 - icu) +iu; = 0, (36) 

where if = ^rju'^^ = rf / n. Note that the parameter H is related to the Coulomb collision 
frequency v^e = rfriQe^ /m and the electron plasma frequency = [jiQe^ jwb)^!'^ as H = 
2(mno//io)(t'ee/'^pe)^- Setting lj = Q — i'y (fi > 0,7 G R), we obtain the dispersion relation 
with respect to the real frequency Q and the damping rate 7, 

f-f + i(l+*' + ^)o^-^ = 0. (37) 



and 



n'(4fi2_C)2-^(C3 + F2) =0, (38) 



where C = 3(P + 1/H) - 1, F = 9[A;2 - l/(2if)] + 1. In this section and AppendixEl we 
use 7 to denote the damping rate. We also have the relation between Q and 7, 

fi2 = 3^2 _ 2^ + p + i_. (39) 

The dispersion relations with various H are shown in Fig. [TJ Note that the determinant 
of the cubic equation (!37|) with respect to 7 is = — (C^ + F'^)/{9 x 36^). If 7 has three 
different real solutions, > 0, i.e., ^^(4^2 _ q2 ^ (i/27)(C3 + F^) < 0, then Q has no 
real solution. Therefore, we have to consider range of the single 7 solution, D.y < 0. This 
range is given by A; > A;crit, where the critical wave number ^crit is defined by i7 = if a 
solution of = exists and fccrit = if there is no solution (see Fig. [T]). These results 
clearly show that the group velocity Vg = dfl/dk is larger than one when H > 3; that is, 
superluminal wave packet propagation is possible (see Appendix [C]). When H > 3.5, there 
are points of = at A; = fc^rit > 0. Figure [2] shows the value of kcrit for each H case. 
The group velocity dVt/dk is infinity dX k ^ kcvit + 0. On the other hand, in the cases of 
H = 1 and H = 2 (see Fig. [U^a) and (b)), the derivative of fl with respect to k increases 
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monotonically, and dfl/dk approaches unity when k becomes infinity (hm^ dQ/dk = 1); 

that is, the gradient remains less than unity as long as k remains finite. Furthermore, we 
prove that dVt/dk < 1 when if < 1.5 (see Appendix [Bll . We find that there is no possibility 
of superluminal propagation of electromagnetic wave in the case of if < 2, while dVt/dk is 
larger than unity in a certain range of k when if > 3. (A detailed investigation produces a 
more strict condition on superluminal propagation of if > 2.3.) 

We now show that matter composed of electrons and positrons with ff > 1 cannot be 
treated as a plasma. This means that superluminal propagation of electromagnetic waves 
is not permitted when the medium is a plasma; and, when it is not, the medium must be 
treated in a different manner. Note that in this paragraph only, we shall use the SI unit 
system. The plasma parameter is given by 



'VP = "A?,= ^(^. (40) 

where T is the temperature of the electron/positron fluids and Ad is the Debye length, 
Ad = ^qT / {ne^) [24]. For a plasma, A^p is (much) larger than unity because charged particles 
are bound to each other when Ap < 1. The frequency of electron-positron Coulomb collisions 
can be written as, 

ne'^lnA ne^ , . 

^ee = 1= = n . /n In A , (41) 



where In A is the Coulomb logarithm and In A' = In A/(6-\/37r) ~ 1/3 24|. Here we used 
InA ~ 10. Then we flnd that H = 2{mn/h){u^Juj^^f < 2(z/ee/ojpe)2 = 2{\nkfne^/{eQTf. 
Finally, we get the relation between ff and Ap, 

HNl < 2(lnA')' ~ ^. (42) 

Therefore, when we consider a plasma (i.e., A^p > 1), we flnd that ff < (lnA')^/A^p < 2/9 < 
1. This clearly shows superluminal propagation of electromagnetic wave is not permitted in 
a true plasma (usually, A'p ^ 1). 

In the above discussion we used the rough approximation InA ~ 10. If In A were greater 
than 100 separately with other variables, if would become larger than 3 and then superlu- 
minal communication would become possible. However, this situation can never be realized 



because there is strict relation between A and Ap as A = GvrA^p (see 



2^, page 24). So, if 



InA becomes larger, then Ap becomes much larger and ff decreases to a value much less 
than unity. 
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FIG. 1: Dispersion relation for electromagnetic waves in the resistive pair plasma for various H. 
The dotted line shows the dispersion relation fcft electromagnetic waves in vacuum. 
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FIG. 2: Dependence of fccrit on H. At k = Aicrit + 0, dO,/dk becomes infinity for each H. 

Next, we discuss briefly the dispersion relation for electromagnetic waves in a uniformly 
magnetized pair plasma. We assume the background plasma is the same as that of the 
previous unmagnetized case except for a uniform magnetic field, B = Bq ^ 0. Using 
the same procedure employed in the previous unmagnetized plasma case, we obtain the 
linearized equations. 

El + -^{Ji X Bo) X Bo = (rj- inu)Ji, (43) 
uj'^Ei + iujJi = k'^Ei, (44) 

where we assume k ■ Vi = and k ■ Ei = to investigate transverse modes. When we 
separate the perturbations of the electric field and current density into two components 
parallel and perpendicular to the background magnetic field Bo, 

El = E\\ + E^, E\\ II Bo, E^ ± Bo, (45) 
Ji = J\\ + J±, J\\ II Bo, J± ± Bo, (46) 

(47) 
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equations (H3ll and (jHj) yield 



-^Jx = (r/-mcj)Jx, (48) 

- uj'^)E^ = iujJ±, (49) 

-^11 = iv " ii^i^)J\\, (50) 

(fc2_cu2)£;|| = iuJ\\. (51) 



(52) 



Finally, we obtain the two dispersion relations, 



(iB^ \ 
r]-iKuj + ] = iuj, (53) 

{k"^ — uj'^){rj — iKoj) = iuj. (54) 

Equations (l53ll and (l5lll are not satisfied simultaneously when i?o 7^ 0. Equation ( IMll is 
the same as that of the unmagnetized pair plasma case. Therefore, we shall investigate the 
dispersion relation ([2SD- 

When we set Co = ku/t], k = nk/r], we have, 

H{cu^ - k'^){iu + - a) =uj^, (55) 

where H = rf / k and a = BlK/ih^rf) = u\/ H {ua = Bq/^/It^ is the Alfven four- velocity). 
Setting ci) = r2 — 27 (r2 > 0, 7 G R), we obtain the dispersion relation for electromagnetic 
waves in a magnetized pair plasma, 

H[{Q^ - 7^ - P) (Q^ -a + -i)+ 2^2^(1 - 27)] = n^- 7^, (56) 
i7[(fi2 _ ^2 _ _ _ 2^(^]2 _ ^2 _ ^ ^ ^ _2^^ (57) 

where we assume i7 7^ 0. Figure [3] shows the dispersion relation for in the case H = 1 
and a = 0.1, 1. The figure clearly shows the group velocity of electromagnetic waves in a 
magnetized pair plasma is less than the light speed in vacuum when H = 1. A detailed 
investigation shows that this is true when < 1 as in the unmagnetized pair plasma case. 
However, when H > 4, the group velocity is larger than the speed of light for some ranges 
of A;. 
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FIG. 3: The dispersion relation of electromagnetic waves in uniform, magnetized pair plasma 
(thick solid lines) and in unmagnetized plasma for comparison (a = 0; thin solid lines), (a) Sub- 
relativistically strong magnetic field case, H = 1, a = 0.1. (b) Relativistically strong magnetic 
field case, H = 1, a = 1. 

IV. CAUSAL RESISTIVE RMHD EQUATIONS 

In the above discussion, we derived the one-fluid equations f|T9l) - fl26l) from the two-fluid 
ones. The one-fluid equations confirmed that the superluminal propagation is forbidden in 
a two-component medium that has a plasma whose plasma parameter greater than unity. 
That is, the one-fluid equations of a pair plasma (fT^ - (|2^ are causal. When we neglect 
the first term of the right hand side in Ohm's law equation (pTI) . which comes from the 
inertial effect of the positron and electron, the term, —iuj/{2v'^^, on the left hand side 
of the dispersion relation (135|) drops out. In this case, the group velocity becomes fg = 
duo/dk = 2A;(4A;2-r/-2)-i/2 > ^hich means the group velocity is greater than the speed of 
light (superluminal). As shown in Appendix [0 when the group velocity is larger than the 
speed of light, superluminal communication would become possible, allowing us to develop 
a device that could send information into the past. However, such a device would destroy 
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the causality of time-ordered events and, therefore, should not be possible. This means 
that, in order to preserve causality, we cannot neglect the inertial term of the electron and 
positron in Ohm's law fl2ip . Recently, several groups performed simulations of resistive 
RMHD including Ohm's law without the electron/positron inertia effect Q, [3]. As shown 
in the above results, unfortunately, all of these calculations are acausal. Here, we propose a 
set of causal resistive RMHD equations in a simple form. 

For simplicity, assuming that \Ah\ ^ h and p+ ~ P-, we obtain the following equations. 
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h 



ne 



V-E 

V B 

OB 



dt 



(UJ + JU) 

pe, 
0, 

-V X E, 



dE ^ ^ ^ 

-^J^V.B, 



0, (58) 
-Vp + peE + J X B{59) 

J ■ E, (60) 

(61) 

(62) 
(63) 
(64) 

(65) 



where 6 = 26{m/e)'^{Q'^ + W) /[p^ — (mQ/e)^] (see Appendix R|) . Here we have to assume 
7' ^ 1 and M = Ui^U" = — 1 in Appendix [A], which means that the relative velocity 
of the electron fluid and positron fluid is nonrelativistic. This condition also preserves 
7 = 1/(1 — V'^y^'^. In a pair plasma, we use 9 = 1. 

A covariant form for these one-component fluid equations (I58|) - fl65|) is as follows: 



9. 



h WU" 



1 



duipU 



0, 



(2ne)2 

U^'Fii - r][r + (f/V,)(l + Q)U^] 



4ne^ 
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n 



(66) 
(67) 

(68) 

(69) 
(70) 
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The difference between equations (l58!)-fl65l) and the RMHD equations used by the pre- 



vious acausal resistive RMHD simulations 
suggested by other articles [l6|,0,Q,ll9|,|2O|. 



21 



5| is mainly in Ohm's law, as expected and 
. The linear analysis of the electromagnetic 



wave in a pair plasma shows that the inertia effect of the electron and positron is essential 
in preserving causality. Here, the electron/positron inertia term of Ohm's law is the right 
hand side of equation (16T!) . If we neglect the change of 7/i/n, Ohm's law simplifies to 



E + V xB--[J--f\pe-V-J){l + e)v] 

7 



d_ 

di 



J + p,V) + V -{VJ + JV) 



(71) 



where k = ~h/{2enf. When H = r]'^/K<l, that is r] < equations dSHD-dSOD, (|62D-(|65D, 
( I7T1) are causal, and thus we call equations (I58l)- (l60l) . (!62l) -( l65l) .( !7Tl) with rj < a/k the "causal 
resistive RMHD" equations. Among the causal resistive RMHD equations, equation ( 17T1) is 
most important; we call it the "causal Ohm's law". When we set B = , equation fl7T|) 
reduces to the simpler Ohm's law. 



E + VxB--[j- -f\p^ - V ■ J)V] 



K 



7 



di 



J + p,V) + V -(VJ + JV) 



(72) 



which is quite similar to the generalized Ohm's law derived by and 21|, but not identical. 



V. EXPECTED PHENOMENA RELATED TO SUPERLUMINAL WAVE 
PACKET 

In this section, we discuss phenomena related to superluminal propagation of electromag- 
netic wave packets that appeared in RMHD simulations that use an acausal Ohm's law with 
H = rf / n > 3. First, we show that it is difficult to detect the superluminal propagation of a 
electromagnetic wave packet in an unmagnetized plasma at rest. For simplicity, we use the 
acausal Ohm's law with k = {H ^ oo). This is just the case of the previous studies with 
resistive RMHD In this case, the dispersion relation becomes that of the telegraphic 

equation, 

uj^ + -uj-k'^ = 0. (73) 
V 

The group velocity of the dispersion relation Vg = dto/dk = k/[k'^ — (2//)"^]^/^ > 1 is 
always greater than the light speed in vacuum. The damping time of the wave is Tdamp = 
l/(— 53(co')) = 21]. The diffusion time of the wave packet is calculated by 



Tdiff 



\D\ \d^Lu/dk^\ 
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where o is the width of the wave packet (see Appendix [Cl equation (IC13P ). The hfe time of 
the wave packet is estimated by 



T 



,T-damp r^isj l + a2(2r7)[P-(2r/)-2]3/2' 
The characteristic propagation length of the wave packet is 



(74) 



^ l + a2(2r7)[P- (2r7)-2]3/2- 



Using N = ka and x = 21] k, we have 



(75) 



Note that the hmitation l/a ^ oo {x ^ oo) means the wave packet propagates to a very 
long distance compared to the scale of the wave packet itself in a highly resistive plasma, 
where the situation is almost the same in a vacuum. However, to detect the superluminal 
propagation of the wave packet, we have to detect the difference between the propagation 
length of the wave packet and that of the light in vacuum, Al = I — t. The difference is 
estimated as 

Al^l/_1_\ N{x' - 1) ^ 22/3iV ^^^^ 

a a\ vj [jY2(^2 _ 1)3/2 + ^2] (^i + ^1 - ^-2^ ~ 3iV4/3 + 22/3- 

Because ^ 1, equation (1751) shows Al ^ a. This means detection of the superluminal 
propagation of the wave packet is difficult in the rest background plasma with a detector of 
ordinary sensitivity. 

When we consider a moving plasma with relativistic speed, propagation of the superlumi- 
nal wave packet changes drastically. Here we consider the wave packet propagating along the 
X direction of a frame {t,x) in a uniform, unmagnetized plasma at rest (see Fig. 111(a)). We 
assume that the wave packet propagates with the group velocity f g > 1 and damps with the 
damping rate 7dmp- Next, we consider a new frame (t',x') moving with velocity vq > l/vg 
relative to the frame {t,x), where the t'-axis and x'-axis in the space-time {t,x) are drawn 
as shown in Fig. Hl^a). The world line of the wave packet is located between the x'-axis and 
X-axis. When we ride on the new frame {t',x'), we see from time inversion arguments that 
the wave packet propagates from the right to the left as shown in Fig. Hl^b). Furthermore, 
the wave packet grows at the rate 7gj.^ = 7dmp \/^~ '^11 ("^^o^g — !)• Here the points A, B, and 
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FIG. 4: Propagation of electromagnetic wave packet in the uniform, unmagnetized plasma, (a) 
Case of rest plasma, (b) Case of relativistic flow of plasma. 

C with respect to the wave packet are identified with those at A', B', and C, respectively. 
This suggests that we have to use the causal RMHD equations (!58l)- (!60l) . (!62l)-(!65l). and 
( I7T|) to avoid such a strange instability of the wave packet, at least in relativistic plasma 
flow, because wave packets propagating in such a flow will grow explosively. Relativistic 
flow exists around the black hole horizon in the Kerr space-time, so artificial radiation of 
electromagnetic wave packets from the horizon will occur in acausal RMHD calculations. 
On the other hand, the same acausal RMHD equations with k, = cause no problem for a 
non-relativistically moving plasma. 

VI. CONCLUDING REMARKS AND DISCUSSION 

We have derived the dispersion relation for electromagnetic waves in a resistive pair 
plasma based on the relativistic two-fluid model, and have shown that the group velocity 
of electromagnetic waves in a resistive plasma is smaller than the light speed within the 
plasma condition (A^p > 1). This shows that superluminal communication is impossible in 
a resistive plasma, thus confirming the causal nature of signals in plasmas. Furthermore, 
the causality condition, H = rf / n < 2(3A^p)^^ < 1, provides an upper limit for electric 
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resistivity in resistive RMHD: 

For simplicity, we assumed that the relative velocity of the positron and electron fluids is 
much smaller than their internal thermal velocities. This is consistent with the assumption 
of linear analysis. In general, however, this assumption is not valid, especially for relativistic 
plasma around black holes. To deal with plasmas where the relative velocity of the electron 
and positron fluids is relativistic, we have to return to the relativistic Vlasov-Boltzmann 
equation with oC.iona. te.n. Q. H Q to obtain the .es,stive tenn ,„ the eausal 
Ohm's law. In such resistivity depends on current density. 

We emphasize that the inertia effect is important for preserving causality, i.e., to forbid 
the superluminal propagation of electromagnetic waves in a resistive plasma. If we neglect 
the inertia term in the generalized Ohm's law (the first term of the right hand side of 
equation (pT!) ). then the resistive RMHD equations give a group velocity of electromagnetic 
waves, fg = duj/dk = 2k/ [Ak'^ — r/"^)^/^, that is greater than the light speed. This shows 
that the inertia effects of electrons and positrons should be considered to preserve causality. 
We therefore proposed a set of causal resistive RMHD equations in section IIVI Numerical 
techniques for simulating "causal resistive RMHD" flow should be developed quickly and be 
applied to astrophysical calculations e.g., energy extraction from a rotating black hole 

rn 

by magnetic reconnection [25]. To perform causal resistive RMHD simulations of a black 
hole magnetosphere, we are to use the general relativistic MHD equations along with the 
causal Ohm's law. 

When we consider matter with a plasma parameter less than one, we cannot use the 
simple two-fluid approximation, because particles in the system are electrically bound to each 
other. Metal, like iron, is an example for such matter. To treat such a relativistic system, 
we ultimately must use relativistic quantum mechanics. However, we have no framework for 
that at present; nevertheless it is an interesting and challenging field for future work. In such 
an unknown framework, the group velocity of electromagnetic waves in any medium should 
not be larger than the light speed (even if the wave damps quickly) to preserve causality. On 
the other hand, within the classical framework where we neglect quantum effects, we would 
show that the group velocity is always equal to or smaller than the speed of light when we 
treat the system properly. Here we cannot use the method of smoothing the electromagnetic 
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field, as is done in traditional particle simulations, because the simple smoothing destroys 
causality. Therefore, numerical calculations may be more difficult than traditional plasma 
particle simulations. It is interesting and important to investigate, in a medium with a 
small plasma parameter, which effects (quantum or classical effects) are more important for 
keeping the group velocity of electromagnetic wave equal to or less than the speed of light. 

In this paper, we considered only a pair plasma; we did not treat an electron-proton 
plasma. However, the similar conclusions also should be drawn for the latter (at lease 
when the plasma is unmagnetized), because the linearized terms of resistive RMHD in 
a pair plasma and in an electron-proton plasma are expected to be similar. It also is 
important to note the differences between RMHD of a pair plasma and in an electron-proton 
plasma. These come from the inequality between the mass ratios of the electron-positron 
and electron-proton. With respect to equations fl58l) - fl60|) . fl62|) - fl65|) . flTTl) . it is expected 
that these equations are similar except for appearance of the second term in the brackets 
on the left hand side of equation ( !59l) . JJ/{2nef', and the term in the brackets on the left 
hand side of equation (17B . 7^(pe — V ■ J)QV . In the electron-proton plasma, the electron 
inertia term with J J is negligible compared to the proton inertia term with UU ^ and 
vanishes because of poor energy exchange between the electron and proton fluids. However, 
in the pair plasma, G is not negligible. Furthermore, the Hall effect disappears in Ohm's 
law fl7T]) in the pair plasma case. Note that all of the terms are nonlinear and that the 
coefficient k of the inertia term in the causal Ohm's law for the electron-proton plasma is 
much smaller than that of the pair plasma (by the ratio of the electron and proton masses). 
It is believed that an accretion disk in a black hole magnetosphere of an AGN will consist 
of an electron-proton plas ma and a corona around the disk and a relativistic jet from AGN 
consist of pair plasma [2611 . Comparison between phenomena in relativistic pair plasmas and 
electron-proton plasmas is both interesting and necessary for understanding the physics of 
black hole magnetospheres where a relativistic jet may be produced. 
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APPENDIX A: DERIVATION OF FRICTIONAL FOUR-FORCE DENSITY 

In this appendix, we derive the friction four-force density between electron and positron 
fluids, whose proper densities are n±. We use and to denote the friction density of 
the electron and positron fluids, respectively. The principle of action-reaction is expressed 
as 

/i + /i = (i = l,2,3), (Al) 

in any inertial frame x'^. When we consider any other inertial frame X'^ = A'^x^ , the 
principle of action and reaction is 

f; + = ^^ui + /-) = ^oUl + /-) = 0. (A2) 
Because 7^ in general, we have /° + /° = 0, or 

n+r-^ 0. (A3) 

Note that + /° = is the law of conservation of energy. We consider the center-of-mass 
frame of the two fluids x*^' — a^x" where the four- velocity of the electron/positron fluids 

satisfies 

n+U'l + n^U'l ^ 0. (A4) 

With respect to the inverse transformation, x^ — b^x"', we have 

n+Ul + n^Ul = Kin+U:;: + n.U'^) - b^M+ + n_7^), (A5) 

where the prime denotes the variable observed in the center-of-mass frame. Using the 
definition — (n+C/^ + ri-U^) / {2n) and 7 = {n+^+ -|- n_7_)/(2n), we have 

K = y (A6) 
In the center-of-mass frame, the spacial components of the friction force density are 

- rl ^ft^ -ma^v,n+n^i^i_ {v% -v^), ( A7) 

19 



where (Tee is the electron/positron colhsional cross section, which is a function of the thermal 
velocity. The average relative velocity of the electrons and positrons, f r, is roughly given by 
the maximum of the thermal velocity and the relative velocity of the two fluids. We write 
the friction four-force density as 

= -ma,,v,{n_^'_n+U^ - n+^'^n_Ut) + 6^/°'. (A8) 

When we use the collision frequency of the electron and positron u^e = c^ee^rn and equation 
[6]), we have 



n ^+7+ + n^j_ 



Next we consider the energy gain rate of the positron fluid in the center-of-mass frame. 
The positron and electron fluids lose the kinetic energy due to friction at the rate, 

= - -^f^n^K, = (l + ^) . (AlO) 



In the above calculation, we employ the principle of action-reaction (ESj), the condition of 
the center-of-mass frame flA4p . and the assumption that the lost energy is thermalized. A 
fraction Q of this thermalized energy (0 < < 1) is distributed to the positron and electron 
fluids, assuming equipartition, and other part is returned to the original fluid. Then the 
energy gain rate of the positron is calculated as 

f% = + (1 - ^)i-K<^) - (All) 

Using the definition of the average four- velocity and four-current density f|T5l) . f|T6l) and the 
friction force density expression (1A9I) . we have 



f+v+i = -m(Teefrn+n_7+7_(f_^ - t;_)f_, 



+i 



2r2+7',_7'e' 



^^[Q^-iyM], (A13) 
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where M = UM", Q = JM", W = J^J". Equations fRT2|) and flXTsD yield 



~(n2M)2-(iLg) 



Here we used 



(A14) 



/ / / tO/ 

7+^+ — = 7 J 



-Q, 



1 

7'^n+7^n_7^ 



12 



-U^'U' = -WU^ = -M, 



u\2 



2e 



The equation of friction four-force density (lA9p reads 

7'e 7 



(A15) 



+ (A16) 



r 



(A17) 



When we introduce the dimensionless factor with respect to the left hand side of equation 



e = ^J^^-^lzm^o. (A18) 



e2 (2mnM)2 - (mQ/e)^ ' 



we finally obtain 



nappVrH 



M ^ ' 



We also calculate the resistive term in Ohm's law 



-r]- 



7 



Q 



J" - —U^(l + Q) 



(A19) 



(A20) 



where 77 = maeeVr/e!^ is resistivity. When we use the collision frequency Vee = Cee'i'r^^, we can 
write 7] = mUee / {ne"^) . 



APPENDIX B: FORBIDDEN RANGE OF SUPERLUMINAL COMMUNICA- 
TION 

We prove that dQ/dk < 1 in the dispersion relation (138!) when H < 1.5. In this appendix, 
we omit the hat over k. The determinant of the cubic equation (1371) with respect to 7 is 



21 



= — (C^ + -F^)/(9 X 36^). If 7 had three different real solutions, it would yield > 0, 
i.e., fi^(4i7^ — C*)^ = {C^ + F^)/27 < would be pure imaginary). Then 7 has only one 
real solution, so that Q has a real solution. When we consider a function of the left hand 
side of equation ([ST]), 



/(7)=7'-7' 



1 + e 



1 

H 



7 



1 

8^' 



we have /(O) = —1/{8H) < and /(1/2) = k'^/8 > 0, and then the single solution of 7 
should beO<7<l/2. Then we have 



because 



37(7 -o) + ^+ 77 = 3 7-77 -- + k' + —. 



From equation (!39l) . we have 



if 



if 



Using equation (1371) . we obtain 

(i7 



dk 



k 



1 - 



372-7 



2 372-27 + i(l + fc2 + l/ii) 



(Bl) 

(B2) 
(B3) 

(B4) 



When H < 3, the denominator in equation (lB4p is positive because of the right side equation 
of equation ( ]B2I) . From equations ( 1B3P and ( lB4p . we have 



dk 2n fi2_3|l__3_+l 



(B5) 



We consider the difference between the numerator and positive denominator of equation 

dESD, 

/ 7^2 1 1 \ / Q ;^2 Q 1 \ 

(B6) 



1 



1 



2 2ii 2 
After some algebraic calculations, we have 



A = k n^--- — + --^-2n -—)- — + - . 



3k' 



AH 



A = ik-n)(n^-k^-l^y^ (4^^ - - :! + 2) + - 7) k-lne + {k-n)^. 

(B7) 



From the left inequality in equation (IBip . we find Q > k when H < 3. And then we have 



(B8) 



^<ik-n)[ — -;-]+{k-n)^-^+[=--j]k~^e. 



2H 



6 
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Using equation (l39l) . we obtain 

A<ik-n)(±- 



6 



k 
2 



H 



5k 



-37^ 



(B9) 



When k > 2/y/2TH, we get ^'^ < P + 1/H < using equation ([BT]). Then we have 

> 5k/2 and A < when H < 3/2. 



On the other hand, when k < 2/V2TH, we obtain Vl^ < k'^ + 1/H < 2h/{2lH) and then 



Vt < 5/V2IH. After some calculations, we have 

'1 



A < (n-k) 



1 kn 3k^\ k 
Q~ 2H^^ T ) ^ 2 



11 



6 A2H 



1 

H 



~4 2~ 



3P _ 3-r^k 
~4 2~ 



< 0, 



(BIO) 



when H < 3/2. Summarizing above calculations, we conclude that A < when H < 3/2. 
This shows that Vg = dVt/dk < 1 when H < 3/2. 



APPENDIX C: PROPAGATION AND DAMPING OF ELECTROMAGNETIC 
WAVE PACKETS 

Here we consider the propagation of a packet of electromagnetic waves in a resistive 
plasma. The wave packet is regarded as an element for communication in the medium. 
First, we use the analytic approximation of a wave packet with a large width. 



1. An analytic approximation solution 

Any variable perturbation of the electromagnetic wave packet in resistive pair plasma, 
/i, is given by, 

/•oo 

= / F{k)e'^''~'^^^^'dk, (CI) 



where F{k) is the Fourier transformation of the variable /i. We take the Gaussian distribu- 
tion of the wave packet to be 



F{k) oc -^e-^ (C2) 



^27r 

where a is the width of the wave packet and fco is the characteristic wave number. The 
initial profile of the variable /i is proportional to exp[— a;^/(2cr^)] ex^i^kox). 
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When l/cr is much smaller than the characteristic scale Ak of the dispersion relation with 
respect to u!{k), we use an approximation 



J ko— J ko— 

where we have expanded to 2nd order in k — ko. When we write u — dk/dLj{ko), D — 
d^u/dk^{ko), 



ikx- 



ko- 



^ifcx-i[w(fco)+|f (fco)(fc-fco)+||^(fco)(fc-fco)' 



f ^tkox-iu}{ko)t 

^ .7 0- 



^-^{a^+iDt)k'^+i{x-ut)k' ^j^l 



ikox—iw(ko)t 



+00 



g- 1(^2 +iDt)k''^ +i{x-ut)k' 



^ikox—iuj{ko)t 



27r 



Here we used a formula, 



TT 



\a'^ + iDt\'^ + ^(a'^ + iDt) 

1 



1 + 



iD*t 



a^ + iDt\^ 



X exp 



^{a'^ + iDt){x-uty 



A - 775 



V2 

X exp 



1 + 



i/coa; — iuj{kQ)t — i 



X exp 



■ 2(a;-Mrt)Mi(o-2_i:);t)t+D,.{(a--nri)2-(mf)2}t 
2{(a2_Dit)2+(I3,t)2} 

2{(a^-At)2+(D,t)2)} 



The width of the wave packet is approximately given by 

Ax - 



(a2 - D,ty + (Ai)^ 



(72 - D;t 



(C4) 
(C5) 



(C6) 



(C7) 



where c is an arbitrary complex constant. When we write u — Ur:-\- iuj and D — Dj. + iD^ 

{ur,Ui, Dr, Di e R), after some algebraic calculations, we have 



(C8) 
(C9) 

(CIO) 



(Cll) 



(C12) 



When t -C „ d.'^d ^. ; the propagation velocity of the wave packet is 

The diffusion time scale of the wave packet, Td, is found from the condition (Aa;)^ = 2a'^, 
which yields 

_2 _2 

(C13) 



T, 



D 



The propagation velocity, t'packet) has meaning only when t <^ or cr ^ ^/t\D\. 
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2. A method of numerical integration and application for imaginary superluminal 
propagation of wave packet 



We show a numerical solution for an electromagnetic wave packet propagation in resistive 
pair plasma, using Simpson's formula. The variable perturbation of the electromagnetic wave 
packet /i given by equations (ICip and flC2l) is calculated as 

3^^(/l)(xe-^°*P(x,^), (C14) 

pko+AI/a 2 

P{x,t) = ^= e-'^(^-^o)'-(^(*^)-^o)*cos(A;x -r](A;)t)dfc, (C15) 

where 70 is the characteristic damping rate of the wave packet. When the number M is 
large enough, the integration becomes the exact value. Usually we set M =4-10, which 
gives precise enough evaluation. To calculate the profile of the wave packet, we evaluate the 
profile function P{x,t). 

Figure shows mathematically the imaginary time evolution of a superluminal electro- 
magnetic wave packet with cr = 50, ko = 0.25 in the pair plasma when if = 4 in equation 
f l36|) . The magnification rate of the variable /i is indicated by the factor beside the ordinate. 
The propagation velocity of the wave packet is found to be f packet ~ 1-4, i.e., superluminal. 
The group velocity Mj. = 1.47 gives a good approximation to the propagation velocity of the 
wave packet. When t = 4, 000 < Tb = 5, 000, the wave packet begins diffuse. The value 
of Td gives a good estimate of wave packet break-down. This numerical calculation clearly 
shows that superluminal propagation of a wave packet is possible if the group velocity of 
electromagnetic wave exceeds the speed of light. 
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FIG. 5: Imaginary time evolution of wave packet superluminal propagation in pair plasma, with 



a = 50, ko = 0.25, H = A. Here /i is normalized by the maximum initial value. 
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